Level statistics of XXZ spin chains with a random magnetic field 
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The level-spacing distribution of a spin-i XXZ chain is numerically studied under random mag- 
netic field. We show explicitly how the level statistics depends on the lattice size L, the anisotropy 
parameter A, and the mean amplitude of the random magnetic field h. In the energy spectrum, 
quantum integrability competes with nonintegrability derived from the randomness, where the XX Z 
interaction is modified by the parameter A. When A ^ 0, the level-spacing distribution mostly 
shows Wigner-like behavior, while when A = 0, Poisson-like behavior appears although the system 
is nonintegrable due to randomness. Poisson-like behavior also appears for A ^ in the large h 
limit. Furthermore, the level-spacing distribution depends on the lattice size L, particularly when 
the random field is weak. 
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I. INTRODUCTION 

Random matrix theories have been successfully applied 
to analysis of the spectra of various physical systems 
such as quantum spin systems strongly corre- 

lated systems^ and disordered quantum systemsi 8 i 9 i 10 i n 
In quantum spin chains, if a given Hamiltonian is inte- 
grable by the Bethe ansatz, the level-spacing distribution 
should be described by the Poisson distribution: 

fpoiO) = exp(-s). (1) 

If it is not integrable, the level-spacing distribution 
should be given by the Wigner distribution: 

ITS f IT S^ \ 

iW*) = yexp [- — )■ ( 2 ) 

In the Anderson model of disordered systems, Pp i(s) and 
-Pwig(s) characterize the localized phase and the metallic 
phase, respectively^ 

The numerical observations should be 
importanti 1 ^ 3 ^! 5 ! 6 ! 7 ! 8 ! 9 ! 10 ! 11 In fact, there has been 
no direct theoretical derivation for the suggested be- 
havior of the level-spacing distribution. Furthermore, 
unexpected behavior has been recently found in the 
level statistics of some XXZ spin chains. 6 Robust 
non- Wigner behavior has been seen in the level-spacing 
distributions of next-nearest-neighbor coupled XXZ 
chains, although they are nonintegrable. The reason why 
it appears is not clear yet, although we have considered 
two possible reasons: extra symmetries or finite-size 
effects. 

In this paper, we discuss the level-spacing distribution 
of a disordered XXZ spin chain so that we may find 
possible clues to the unexpected behavior of the XXZ 
spin chains. We consider specifically the spin-i XXZ 
spin chain with random magnetic field, where quantum 



integrability competes with nonintegrability due to the 
randomness. The Hamiltonian on L sites is given by 

L L 

3=1 i=i 

(3) 

where S a — (l/2)a a and (a x , a y , a z ) are the Pauli ma- 
trices; hj is random magnetic field along the z axis at 
site j; the periodic boundary conditions are imposed. 
The random magnetic field hj J s are uncorrelated ran- 
dom numbers with a Gaussian distribution: (hj) = 
and (h n h m ) — h S nm . Here we recall that the system 
is integrable when h = 0, while it is not when h ^ 0. 

We show explicitly how the level-spacing distribution 
P(s) depends on the lattice size L, the anisotropy param- 
eter A, and the mean amplitude of the random magnetic 
field h. For some special cases of the Hamiltonian ©, the 
level statistics has been discussed. The present numeri- 
cal results should make explicit connections among them 
and even extend them. When A^0, -P(s) mostly shows 
Wigner-like behavior. For the Heisenberg case (A = 1), 
it has been shown that P(s) coincides with Pwig(s)- 11 
When A = 0, Poisson-like behavior appears although the 
system is nonintegrable due to the randomness. However, 
the result is consistent with the Anderson localization of 
one-dimensional (ID) systems. The Wigner-like behav- 
ior for A ^ suggests that the Anderson localization in 
ID systems should be broken by the interaction such as 
the XXZ coupling. Here we note that an analogous phe- 
nomenon has been observed in the 2D Anderson model 
with electron interaction^ In the large h and small h 
limits, Poisson-like behavior appears again for A / 0, 
which is consistent with the numerical results for spin- 
glass clusters*^ the open-boundary Heisenberg chain^ 
and the 3D Anderson models Furthermore, we find that 
P(s) depends on the lattice size L, particularly when the 
random field is weak. 
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TABLE I: Matrix size of the largest subspace and the number 
of samples calculated in this work for each lattice size. 



Lattice size 


Matrix size 


Number of samples 


8 


70 x 70 


10000 


10 


252 x 252 


3000 


12 


924 x 924 


1000 


14 


3432 x 3432 


900 



There is another motivation for the present study. The 
symmetry of the integrable quantum XXZ spin chain 
can be nontrivial. An extraordinary symmetry appears 
for special values of A: the XXZ Hamiltonian commutes 
with the sl2 loop algebra when q is a root of unity, where 
q is defined by A = (q + l/q)/2H The loop algebra is 
an infinite-dimensional Lie algebra, and the dimensions 
of degenerate eigenspaces are given by some exponen- 
tial functions of the system size, which can be extremely 
largeii2*i£ It should therefore be nontrivial how the large 
degeneracies are resolved by nonintegrability. 



II. NUMERICAL PROCEDURE 

In the Hamiltonian ©, total S z is conserved. The 
eigenstates with different S z are uncorrelated. Therefore, 
we consider only the largest subspace S z — 0. The largest 
sectors for the lattice size L =8, 10, 12, 14 have 70, 252, 
924, 3432 eigenvalues, respectively (see Table |IJ. 

To find universal statistical properties of the Hamilto- 
nians, one has to deal with unfolded eigenvalues instead 
of raw eigenvalues. The unfolded eigenvalues are renor- 
malized values, whose local density of states is equal to 
unity everywhere in the spectrum. In this paper, the 
unfolded eigenvalues Xi are obtained from the raw eigen- 
values Ei in the following method. Let us define the 
integrated density of states as 



(4) 



Here 9(E) is the step function and N is the number of 
the eigenvalues. We choose some points of coordinates: 
(Ei,n(Ei)) for i = 1, 21, 41, • • • , N. The average of inte- 
grated density of states (n(E)) is approximated by the 
spline interpolation through the chosen points. The un- 
folded eigenvalues are defined as 



Xi = (n(Ei)). 



(5) 



The level-spacing distributions are given by the proba- 
bility function P(s), where s — Xi+i — X{. 

We have calculated 10000, 3000, 1000, 900 samples of 
P(s) for L =8, 10, 12, 14, respectively (see Table 0, 
and averaged the samples for each L. To calculate the 
eigenvalues, we have used standard numerical methods, 
which are contained in the LAPACK library. 




FIG. 1: Level-spacing distributions for L = 14, h/J = 0.5, 
A =0, 0.1, 0.2, 0.5, 1. Broken line, the Poisson distribution; 
solid line, the Wigner distribution. 



III. LEVEL-SPACING DISTRIBUTION 

Depending on the anisotropic parameter A (0 < A < 
1), the level-spacing distribution P(s) changes between 
the Wigner distribution Pwi g (s) and the Poisson distri- 
bution Ppoi(s) as shown in Fig. ^ where L — 14 and 
h/J = 0.5. When A = 0, P(s) almost coincides with 
Ppoi(s) although the system is nonintegrable due to the 
random magnetic field. As A increases, P(s) rapidly 
changes to P\vig(s)- 

Let us explain the Poisson-like behavior of A = 
in terms of the Anderson localization. The Hamilto- 
nian @ can be mapped into a model of interacting ID 
free fermions under random potential: 



n 



j 



■ 2 



L-l 



5>}<*+i + c] +lCj ) (-l) M (c[ Cl + c\c L ) 



4=1 



A 7 I> c 4 c 4 c 4+i c 4+i ~ 4c 4 c 4 + 1) 



4=1 



(6) 



Here, L is the number of sites; M is the number of 
fermions; cj and Cj are the creation and annihilation op- 
erators of fermions on the jth site, respectively. And the 
Anderson model of noninteracting disordered fermions is 
given by 



U = z2 £ 3 c ] c J 



c]Ci), 



(7) 



where Ej is the random potential at the jth site; V is a 
constant hopping integral; denotes summation over 
nearest-neighbor sites. One can find that Eq. |JBJ for A = 
corresponds to Eq.Q of the ID case. It is known that 
localization always occurs in the ID case, while the 3D 
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Anderson model has the metallic phase and the localized 
phase. Here we recall that the metallic phase corresponds 
to Pwig(s) and the localized phase to Pp i(s). Thus, the 
observed Poisson-like behavior for A = is consistent 
with the Anderson localization. 

Let us discuss a consequence of the Wigner-like behav- 
ior for A / 0. The Hamiltonian namely Eq. 
for A / corresponds to interacting ID fermions un- 
der random potential. Thus, the Wigner-like behavior 
of P(s) for A ^ might suggest that the interaction 
among fermions should break the Anderson localization 
in ID chains. The suggestion could be consistent with 
the observation on the 2D interacting lattice fermions 
where the level-spacing distribution changes from Pp i(s) 
to fWig(s) as the electron-electron interaction increases 
from zero^ 

We now discuss how the level-spacing distribution P(s) 
depends on the random magnetic field h. We consider 
only the case of A / 0. In Fig.[5]the graphs of P(s) are 
shown for some values of h/ J and L, where A = 0.5. We 
first consider the case of large h. As h/ J increases from 
the value of 0.5, we observe that the form of P(s) changes 
from Pwig( s ) to Ppoi(s). The observation suggests that 
the effect of random magnetic field on each site should 
become larger than that of the correlation between adja- 
cent spins, as the random field h/J increases. The spins 
should become more independent of each other as h/J 
increases, since the effect of correlation decreases effec- 
tively. Thus, the Poisson-like behavior of P(s) should ap- 
pear in the limit of large h/J. Similar shifts from Pwig(s) 
to Ppoi(s) as randomness increases have been discussed 
for the 3D Anderson model^ the spin-glass clusters*^ 
and the open-boundary Heisenberg chain^ 

For the case of small h, the level-spacing distribution 
P(s) strongly depends on the system size L, and the be- 
havior of P(s) is dominated by finite-size effects. In Fig. 
2, we observe that the form of P(s) for h/J = 0.1 is 
different from that of the standard Wigner distribution, 
particularly when L is small. When L is small, random 
magnetic field is irrelevant to energy levels if it is smaller 
than the order of 1/L. In fact, energy differences should 
be at least in the order of 1/L, and random magnetic field 
can be neglected if it is much smaller than some multiple 
of 1/L. Thus, for the case of small h, the level statistics 
should show such a behavior as that of h = 0. In fact, the 
Hamiltonian for h = is the integrable XXZ spin chain, 
which should have Poisson-like behavior. Furthermore, 
the integrable XXZ Hamiltonian at A = 0.5 has the sh 



loop algebra symmetry^ and the level-spacing distribu- 
tion should show a peak at s = 0. 6 In Fig. 2, the graph 
of P(s) for L = 8 and h/J= 0.1 suggests such behavior. 

Let us discuss finite-size effects on the level-spacing 
distributions. In order to observe the size dependence 
of P(s) clearly, we employ the following parameter, rj — 
J° [P(s) - P wig (s)]ds/ f*°[P Poi (s) - Pwig(s)]^, where 
So = 0.4729 • • • is the intersection point of Pp i(s) and 
Pwig(s)' 10 ' 12 Thus, we have r\ = when P(s) coin- 
cides with Pwig(s), and r\ = 1 when P(s) coincides with 
Ppoi(s)- In Fig. 3, the value of r\ for h/J — 0.1 strongly 
depends on the lattice size L. Moreover, we observe that 
as L increases, r\ decreases for h/J — 0.5, while r\ in- 
creases for h/ J = 2. The observation suggests that 77 ap- 
proaches either the value or 1 as L increases. In other 
words, it should become more definite whether P(s) has 
Wigner-like behavior or not, as the system size becomes 
large. 



IV. CONCLUSIONS 

In conclusion, we have calculated the level-spacing dis- 
tributions of finite spin-i XXZ chains under random 
magnetic field, and shown how the level-spacing distri- 
butions change between the Poisson distribution and the 
Wigner distribution depending on the lattice size L, the 
anisotropy parameter A, and the mean amplitude of the 
random magnetic field h. For A = 0, the level-spacing 
distribution P(s) almost coincides with the Poisson dis- 
tribution although the system has the randomness. As A 
increases from zero, P(s) rapidly shifts to the Wigner dis- 
tribution. The behaviors of P(s) have been explained in 
terms of Anderson localization. For A ^ 0, P(s) strongly 
depends on L when h is small. When L is finite, P(s) 
should show Poisson-like behavior in the small h limit. 
In the large h limit, however, P(s) should become close 
to the Poisson distribution independent of L. 
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FIG. 2: Level-spacing distributions for L = 8, 10, 12, 14, /i/J =0.1, 0.5, 1, 2, A = 0.5. Broken lines, the Poisson distribution; 
solid lines, the Wigner distribution. 
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FIG. 3: Dependence of the parameter 77 on the lattice size L 
for h/J =0.1, 0.5, 1, 2 and A = 0.5. r\ = corresponds to the 
Wigner distribution and 77 = 1 to the Poisson distribution. 
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